Spin-dependent quasiparticle reflection and bound states at interfaces with itinerant 

antiferromagnets 
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We present a formulation of the quasiclassical theory of junctions between itinerant antiferromag- 
nets (AF) and s-wave (sSC) and d-wave superconductors (dSC). For the simplest two-sublattice 
antiferromagnet on a bipartite lattice, we derive Andreev-type equations and show that their so- 
lutions lead to a novel channel of quasiparticle reflection. In particular, quasiparticles in a normal 
metal with energies less than or comparable to the antiferromagnetic gap experience spin-dependent 
retroreflection at antiferromagnet-normal metal (AF/N) transparent (100) and (110) interfaces. A 
relative phase difference of vr between up spin and down spin quasiparticle reflection amplitudes is 
shown to lead to zero-energy interface bound states on AF/sSC interfaces. For an sSC/AF/sSC 
junction, these bound states are found to be split, due to a finite width of the AF interlayer, and 
carry the supercurrent. At AF/dSC interfaces we find no zero-energy bound states for both interface 
orientations we considered, in contrast with the case of (110) impenetrable surface of a dSC. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



Introduction. The study of electronic properties of 
superconductor-magnetic interfaces has enjoyed a renais- 
sance in recent years, with increased interest in techno- 
logical applications and advances in theory. Ferromag- 
netic layers can spin polarize quasiparticle currents and 
Zeeman split surface densities of states, with possible ap- 
plications in spintronics. SFS junctions have been shown 
to display — tt transitions with varying temperature, 
width, or orientational structure of magnetization of the 
ferromagnetic interlayeri2iMiSiSiL2iS. The quasiclassical 
theory of superconductivity, complemented with appro- 
priate boundary condition s--^-?'— , has proven very suc- 
cessful in dealing with problems like superconducting- 
ferromagnetic interfaces. No quasiclassical theory of the 
antiferromagnetic state and, in particular, AF/SC inter- 
face exists, however. 

On the other hand, there are many situations 
of fundamental and practical interest which involve 
antiferromagnet-superconductor interfaces (AF/SC). Ini- 
tial results of experimental investigations of proximity 
and Josephson effects through this type of interfaces have 
been obtained only recentlj^. Some unusual properties 
for specific antiferromagnetic Josephson weak links have 
been studied theoretically for barriers built on a giant 
magnetoresistance multilayer or made of doped mangan- 
ites in the metallic A-phasei^. Many of the properties of 
high-temperature superconducting (HTS) cuprate mate- 
rials are also thought to result from a competition be- 
tween antiferromagnetic and superconducting order, and 
there are many naturally occurring situations and possi- 
ble devices which might involve such boundaries. These 
include interfaces of insulating and highly doped cuprates 
or SC/AF/SC junctions, HTS grain boundaries where 
antiferromagnetism may play a role as a surface state, 
and the antiferromagnetism which has been observed in 
HTS vortex cores. 

Quasiclassical equations describe spatial variations of 



physical quantities taking place over length scales which 
are long compared to atomic distances. The boundary 
conditions match slowly varying quantities at interfaces, 
where the equations themselves do not apply. We ar- 
gue below that the assumption of slow variation of the 
sublattice magnetization allows a formulation of the qua- 
siclassical approach to itinerant antiferromagnets, which 
we study with mean field theory. We have found a novel 
spin-dependent channel of normal metal (N) quasipar- 
ticle reflections from AF/N interfaces, associated with 
the structure of the antiferromagnetic order parameter 
as well as with a magnetic periodicity along the interface 
with component Qy of the antiferromagnetic wave vector. 
Parallel to the interface, components of crystal momenta 
for outgoing and incoming quasiparticles differ by Qy in 
the novel reflection channel, while differences between 
components normal to the interface depend on the shape 
of the Fermi surface. Specific channels of quasiparticle 
reflection have been theoretically found also in charge 
density waves (CDW) - normal metal junctionsi^ii^iiSiii. 
An excess quasiparticle reflection below the CDW gap 
in NbSca has been experimentally identified in Ref. Il8l 
We show below, for the simplest cases of (100) and (110) 
interfaces, that the spin dependence of quasiparticle re- 
flections from boundaries with itinerant antiferromagnets 
results in dramatic consequences which do not appear at 
the nonmagnetic interfaces mentioned above, however. 
An important distinctive feature of the quasiparticle re- 
flection we study below is the phase difference tt of reflec- 
tion amplitudes taken for spin-up and spin-down quasi- 
particles in the normal metal at an AF/N interface. As 
a result of this "7r-shift", zero-energy interface bound 
states arise on AF/sSC transparent interfaces. By con- 
trast, a d-wave superconductor, which in the case of a 
free (110) surface manifests Andreev zero energy bound 
states, shows no zero-energy states at AF/dSC transpar- 
ent interfaces. For an sSC/ AF/sSC junction, the zero- 
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energy states are found to be split and carry a supercur- 
rent. They can form fow-temperature anomalies in the 
Josephson critical current through sSC/AF/sSC junc- 
tions. They also give rise to the zero-bias anomaly of 
the conductance of N/AF/sSC junctions. 

Quasiparticle reflection at a smooth AF interface. We 
begin by considering electrons hopping on a square lat- 
tice with basis vectors a and b (lattice constant a), with 
superconducting pairing A*-' and magnetization m*. 

+ J2m'ini^ -n,i)- (1) 

i i 

We assume nearest neighbour hopping, and consider 
either s-wave pairing A'-' = ~'Vs{ciiCi'^)Sij — SijAl 
or (i-wave pairing A*-'' = -Vd{ciiCjf) = A''^S\i^j\,i 

such that A^'=^° = -A^*^^ The magnetization to' = 
~(Kn/2)(n.i| — riii) as a function of site i is assumed to 
exhibit antiferromagnetic order, such that for the infi- 
nite lattice in the absence of local perturbations, = 
(^—lya+jb^^ In the more general case, the sublattice 
magnetization to is assumed to be slowly varying function 
on scales of the lattice spacing a. We assume always that 
TO is nonzero only on one semi-infinite half-space, while 
A may be nonzero on the other; thus there is generally 
no additional potential barrier between the two systems, 
although we discuss the possibility. 

The normal state electron band ^(k) = — /i— 2i(cos fca + 
cos kb) and the respective Brillouin zone is spanned by 
ka,b G [— 7r,7r], where momenta are given in units of 
a~^. In describing quasiparticle reflections, it is con- 
venient also to work in a coordinate system where x 
and y describe coordinates perpendicular and paral- 
lel to the interface, respectively. For a (100) inter- 
face, the band and the zone have the same forms in 
the y-coordinates. For a (110) interface, however, 
we have ^(k) = —fi — 4i cos(fc^/'\/2) cos(fcy/A/2) and 

kx € [— a/Stt, v^tt], ky e [— 7r/V2, 7r/-\/2], on account of 
the periodic conditions along the surface. 

A specific feature of the antiferromagnetic state, which 
is important in the derivation of quasiclassical equa- 
tions, is the presence of rapidly oscillating term to-' = 
(^—lya+jbjji — exp{iQj)m, which results in slowly vary- 
ing Andreev amlitudes with wave vector fe -I- Q in the 
equations for the amplitudes with the wave vector k. In 
the case in question, 2Q coincides with a basis vector 
of the reciprocal lattice of the nonmagnetic crystal. For 
this reason the quasiclassical equations can be written 
for pairs of entangled quasiparticle trajectories k and 
k + Q. As usual in the quasiclassical theory of super- 
conductivity, we require that = ^Vp/As^d ^ a, and 
similarly require that the magnetic "coherence length" 
= hvF/\m\ 3> a. We also assume that the devia- 
tion from half-filling in the antiferromagnet is not large 

<C ej?. Otherwise the antiferromagnetic state would 



be unstable within the framework of a generic Hubbard- 
like model. Then fi should be included directly in the 
quasiclassical equations, not in the rapidly oscillating ex- 
ponentials. Under this condition, the nesting relation is 
valid with quasiclassical accuracy and energies of normal 
state quasiparticles with momenta kp and kp + Q both 
lie on the Fermi surface. 

It is now convenient to collect into a Nambu 4- 
spinor the Andreev amplitudes il^J = {uja{kF),Uja{kp + 
Q),Vjg{kp),Vja{kp + Qj), and to define the Pauli ma- 
trices Pa, Ta in space of two quasiparticle trajectories (fc, 
k + Q) and in particle-hole space respectively. 

Then the Andreev equations take the form: 

^-/iTspo - iT3P3Vp,x-^ + crm{x)Topi + A(x)^ ■0(a;) = 

= £^(x).(2) 

Here vp^ — ^^r^\ „ is the Fermi veloc- 

-^'^ ak^ \fj,=Q 

ity at half-fiUing, A{x) — As{x) + Ad{kp,x), 
As{x) = apQAs{x)— + o-poA*(a;) — , Ad{kp,x) = 

aAd{kp,x)p3— + <7A*i{kp,x)p3 — . Quasiclassical 

Eqs.lO do not apply in vicinity of quasiparticle mo- 
menta where vp,x — 0. In particular, they do not 
apply near saddle points of quasiparticle energies where 
Van Hove singularities of the normal metal density of 
states take place. Since we will be interested mostly 
in transport across the interface, where the additional 
factor vp_x arises, these momenta do not contribute to 
the results noticeably and the conditions turn out not 
to be restrictive. Eqs.(|2Jl are written for a quasiparticle 
spin direction a (cr = —a = ±1), with the quantization 
axis along the magnetization at site ja ~ jb = 0. The 
quasiclassical equations for the superconducting side can 
be formulated also for n ep. The parameter fi does 
not enter quasiclassical equations in this case and leads 
to a large Fermi velocity mismatch on the boundary. 
We have presented only the result |(5J) valid for /i <C sp, 
since we discuss primarily the case of a small Fermi 
velocity mismatch between a superconductor and an an- 
tiferromagnet. Eq.lPJ also implies that no subdominant 
pairing channels are important in the problem, which 
we solve below with step-like (nonselfconsistent) profiles 
of the order parameters. 

Normal/ AF interface. We now apply Eqs.® to the 
problem of electron refiection at AF/N interfaces. The 
initial and final states of an electron reflected from AF/N 
interface belong to the bulk of the normal metal. We 
find that the incoming electron with a momentum kp 
acquires in a reflection event the momentum kp + Q. 
The antiferromagnetic wave vector on the square lattice 
is Q = (tt, tt), with respect to the crystal axes. In the 
X, j/-coordinate system, Q = (tt, tt) for (100) interface and 
Q = (\/27r,0) in the (110) case. The quasiparticle veloc- 
ity Va.b — 2t sin fca.b chaugcs its sign when the momentum 
varies by Q. Hence, quasiparticles experience retrorefiec- 
tion at AF/N interfaces. Since vp,y = for a (110) inter- 
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face at half-filling, in this particular case quasiparticles 
move along the interface normal, so that retroreflection 
and specular reflection coincide. We note that, beyond 
the quasiclassical approximation, Q is not a difference 
between incoming and outgoing momenta anymore and 
the condition for retroreflection is broken. 

For energies below the quasiparticle gap in the antifer- 
romagnet we find the following reflection amplitudes for 
electrons at the AF/N interface (for both orientations): 

rAF,a = , \^l + £\<m. (3) 

am 

Reflection amplitudes for quasiparticle energies above the 
gap + £| > m are obtained from Eq.© with the sub- 
stitution i^Jm^ — (m + ~^ sgn{ii + e)\/(Ji + e^^^^rn? . 

The retroreflection at an AF/N interface is analogous, 
to some extent, to Andreev reflection at an SC/N inter- 
face. In particular, bound states arise in AF/N/AF meso- 
scopic systems for energies below the antiferromagnetic 
gap, analogously to Andreev subgap states in SC/N/SC 
systems. The eigenenergies of these states are determined 
by the equation e + /i = ^-^^((ys + 7m), with integer n 
and ip — arccos[(e + //)/to]. For sufficiently large width 
d 3> of the normal metal layer, and considering quasi- 
particles with subgap energies |£ + /^| ^ m, we obtain the 

simple, explicit expression Sn + fJ- = {n + As 

is known, in SC/N/SC systems electrons and holes with 
the same momenta and opposite velocities form Andreev 
bound states, which carry electric current. However, due 
to different nature of AF and SC order parameters, in 
AF/N/AF systems electrons with momenta k and k + Q 
and opposite velocities are coherently entangled and the 
bound states themselves do not carry electric current. 

Until now we have not discussed nonmagnetic channels 
for quasiparticle reflection, which are present both on ac- 
count of Fermi velocity mismatch in a normal metal and 
an antiferromagnet and/or due to potential barriers at 
the interface. Potential barriers themselves form specu- 
lar quasiparticle reflection, as usual. If the characteristic 
parameters t and fi in two identically oriented half-spaces 
are close to each other, the reflection coefficient of con- 
ventional specular reflection arising due to a mismatch of 
Fermi velocities is of order (m/i)^. Furthermore, crystal 
momentum can change in a reflection process by a recip- 
rocal crystal vector along the surface. Due to a difference 
between reciprocal crystal vectors at the interface and in 
the bulk, specific crystal periodicity along a particularly 
oriented surface or interface can result in additional chan- 
nels for quasiparticle reflection^^. 

AF/SC interface. Eqs.^ apply also to AF/SC inter- 
faces. The solution of these equations in the absence 
of potential barriers and/or a Fermi velocity mismatch 
shows zero-energy interface states in the system. These 
states take place for arbitrary relation between m ^ ep 
and As ^ Ei? and for < m, which garantees the exis- 
tence of the antiferromagnetic gap for electrons and holes. 
The zero-energy states at the AF/sSC interface arise as a 



combined effect of Andreev reflection from the supercon- 
ducting halfspace and the antiferromagnetic retroreflec- 
tion from the antiferromagnetic side. The origin of the 
zero-energy surface states is closely connected with the 
magnetic properties of the quasiparticle reflection. This 
can be easily seen in the particular case m 3> , where 
quasiparticle reflection from the antiferromagnet can be 
described with reflection amplitudes l|2Jl for the AF/N 
interface. According to Eq.ljJl, antiferromagnetic order- 
ing results in opposite signs of reflection ampltudes for 
electrons with spin up and down. Quasiparticles with en- 
ergies below the antiferromagnetic gap do not penetrate 
in the bulk of the itinerant AF, rAF,a = exp(i0o.). This 
takes place in the absence of potential barriers, i.e. for 
a transparent interface. Quasiparticles in the supercon- 
ducting halfspace can be described in terms of standard 
Andreev equations for Andreev amplitudes iff [kp) = 
(ttif (fep), WixC^F), "^iT (^f), "Cjj. (fc_F)) complemented with 
the boundary conditions ^°%{ky + Qy) = Sipl"Q{ky). 
The Andreev amplitudes V'™ contain solutions for quasi- 
particles moving towards the interface {vf,x < 0), in 
contrast with i/;""*. The S'-matrix for the AF/N bound- 
ary takes the form 5* = (^t'afj — ^ ^ + ^af.i — 2~~^)^°' 
Thus, the problem of a superconducting halfspace with 
AF/SC interface can be formulated for quasiparticles be- 
low the antiferromagnetic gap in a form identical to that 
obtained for an impenetrable ferromagnetic surface. The 
spectrum of Andreev bound states es = ±Acos(6/2), 
where = 0| — is well studied for the latter 
probleroii. As follows from Eq. © , for the antiferromag- 
netic boundary = tt. This leads, indeed, to zero-energy 
surface bound states on the AF/sSC interface. 

Let now a d-wave superconductor make a junction with 
an antiferromagnet without potential barriers and/or 
substantial Fermi velocity mismatch. The AF/dSC in- 
terface is closed for quasiparticles in the superconduc- 
tor with energies below the antiferromagnetic gap. How- 
ever, no low-energy surface bound states with \sb\ -^m 
(in particular, no zero-energy states) form in the case in 
question for any interface-to-crystal orientation. This fol- 
lows directly from solutions of Eqs.Q and the continuity 
of Andreev amplitudes at the boundary. This can be also 
explained qualitatively as a combined result of the anti- 
ferromagnetic retroreflection and the change of sign of the 
d-wave order parameter. The difference tt between phases 
of reflection amplitudes © can be effectively ascribed to 
the variation of the phase of the order parameter also. 
In order to see this, one can introduce auxiliary quanti- 
ties u™(x, fcF,£) = ^^(a;, fef, £)e^*'^®/^, w^"(x, A;_f, e) = 
Vs{x,kF,e)e^'^'^/'^ into Andreev equations and boundary 
conditions, taken for the outgoing momentum kp- An- 
dreev amplitudes for incoming momentum kp are kept 
unchanged. Then the problem becomes, formally, iden- 
tical to the one for magnetically inactive impenetrable 
boundary and the effective order parameter for the outgo- 
ing momenta /^ef },a{xfkp) = e~'^'^^A{x, kp). Since the 
d-wave order parameter Ad{kp,Xi) — 2 A]^^"' {cos ka — 
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cos kb) changes its sign when the wavevector changes by 
Q, we obtain Ad{x, kF) = —'Ad{x,kp) for kp = kp + Q 
Thus, an outgoing quasiparticle sees an effective super- 
conducting order parameter with an additional phase 
TT — as compared with the order parameter for the 
incoming trajectory. For the antiferromagnetic interface 
Q = TT and the total phase variation of the effective order 
parameter in a reflection event vanishes. 

Interface potential barriers open a channel of specular 
reflection. A combined description of specular reflection 
and retroreflection includes more complicated boundary 
conditions, where S'-matrix connects V'o"*(^y): '0o"*(^y + 
Qy) with ij'a" {ky) , iPq" {ky + Qy) and contains both off- 
diagonal and diagonal components in momentum space. 
This modifies the effects of retroreflection. The zero- 
energy surface states will be split at an AF/I/sSC inter- 
face and reach the value of the superconducting gap in the 
limit of impenetrable insulating interlayer. The opening 
of conventional channels of reflection on the AF/I/dSC 
interface will result in interface subgap states. For a (110) 
interface orientation, the subgap states will evolve with 
decreasing barrier transparency to the well known zero- 
energy surface states in d-wave superconductors. For 
large potential barriers, the antiferromagnet weakly splits 
the zero-energy states. A mismatch of Fermi velocities 
similarly results in conventional channels of reflection. 
For a small mismatch, the energies of subgap states on 
the AF/dSC interface will be situated close to the super- 
conducting gap. If the mismatch (e.g. the magnetiza- 
tion) increases, the subgap states on the (110) interface 
will move towards lower energies. We believe this to be 
the origin of the subgap states on the AF/dSC (110) 
interface obtained recently in Ref. on the basis of nu- 
merical studies for a large mismatch of Fermi surfaces. 

Josephson current through SC/AF/SC junction. The 
zero-energy surface states arising on the AF/sSC inter- 
face will be split for an AF interlayer with finite width 
I. The split can be considered as an effect of tunneling 
between zero-energy states on two boundaries of the AF 
layer. If no potential barriers are present on the bound- 
aries and ? <C 'Csj "T- ^ As, we find the following ener- 
gies for interface states: Sb = ±\/^|As cos(x/2)|, where 
D{ky) — AK{ky){K{ky) + 1)~^ is the transparency of 
the N/AF/N junction and K{ky) = eyip{2ml/\vF,x{ky)\). 
These states carry the Josephson current 

Tdfc _ y V:D|A,|cos^ 
J= j ^2ev^|A,|sin|tanh — 2. ^ (4) 

-Tr/2 



which differs from the Ambegaokar-Baratoff result. In 
the particular case of large interlayer width, K,D <^ 1, 
there are low-energy states in the junction which result in 
low-temperature anomalous behavior of the critical cur- 
rent. This behavior is similar to what can happen in 
tunnel junctions with d-wave superconductors or S/F/S 
junctions with low-energy interface states2*ii*2i. 

Conclusions. We have developed the quasiclassical 
theory of itinerant antiferromagnets under the assump- 
tion of slow spatial variations of the sublattice magne- 
tization and near-nesting of the metallic band. The 
quasiclassical equations, written for pairs of entangled 
quasiparticle trajectories, provide a powerful tool for 
study of a host of new problems related to interfaces 
with antiferromagnetic materials. As important appli- 
cations of this formalism, we have shown that quasipar- 
ticles in normal metals are retroreflected from the anti- 
ferromagnet, and studied bound state formation at in- 
terfaces with superconductors. We have further shown 
that Andreev zero-energy surface states are formed at 
an interface between an antiferromagnet and an s-wave 
superconductor. These bound states do not arise on in- 
terfaces with a d-wave superconductor for any interface 
orientation. We demonstrated that the Josephson cur- 
rent in sSC/AF/sSC junction is strongly influenced by 
the low-energy interface states and differs from the con- 
ventional Ambegaokar-Baratoff result. In analogy with 
other situations'^, one can expect that, after introduc- 
ing some effective order parameters, the self-consistent 
calculations do not modify low-energy bound states and, 
in particular, their contribution to the Josephson cur- 
rent. Our formalism in the present paper is based on the 
Andreev equations, but we have derived also the Eilen- 
berger equations for itinerant antiferromagnets and re- 
spective boundary conditions. Based on these equations, 
we plan to explore with self-consistent calculations prox- 
imity effects for AF/SC interfaces as well as the problems 
of the Josephson current and the low- bias conductance in 
a longer work. 
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